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. We show how, contrary to physical intuition, thermal fluctuations of physical 

states having a non- vanishing baryonic number can be fully neglected in the 
thermodynamics of QCD at any physical temperature. We also discuss on 
o the consistency between our results and the CPT theorem. The last part of 

this article is devoted to comment some interesting physical features which 
follow from this result. 
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The partition function of a quantum field theory at zero temperature is 
dominated by the contribution of the vacuum state, i.e., the eigenstate of the 
quantum hamiltonian having minimum energy. By increasing the physical 
temperature, contributions to the free energy density of excited states become 
more and more relevant and in the infinite temperature limit the system 
fluctuates randomly and all the states have the same probability, no matter 
what energy they have. 

In QCD, the gauge theory describing the strong interaction of particles, 
the SU(3) gauge symmetry of the model implies that physical states are 
made up of a multiple of three number of quarks. Hence we have baryons 
and mesons in the spectrum of this model and a natural expectation is that 
thermal fluctuations of both kind of particles will be relevant at finite tem- 
perature. Notwithstanding that we will show in this letter how thermal fluc- 
tuations of physical states having a non-vanishing baryonic charge can be 
fully neglected in the thermodynamics of QCD at any physical temperature 
[p]]. The last part of this article will be devoted to discuss some interesting 
physical features which follow from this result. 

1. Analytical results 

Due to a global U(l) symmetry of QCD, the vector current is a 

conserved current. The integrated fourth component of this current 



defines a dynamical invariant, the baryonic charge or quark number operator. 
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(1) 



After quantization the quark number operator commutes with the hamilto- 
nian and physical states can therefore be characterized by its energy and 
baryonic charge. Hence the Hilbert space of physical states can be decom- 
posed as a direct sum of Hilbert subspaces, each one of them corresponding 
to a given eigenvalue of the quark number operator. 

Taking into account this decomposition of the Hilbert space of physical 
states we can write the partition function of QCD at finite temperature 

Z = Tr (e"T ) (2) 
as a sum of canonical partition functions 

^ = £Tr fe (e-f), (3) 
k 

each one of the Tr^ ( e ^) being the partition function computed in the 
baryonic sector of the Hilbert space corresponding to a given eigenvalue k of 
the quark number operator. 

In order to establish the connection between this hamiltonian approach 
and the path integral formulation in the lattice regularization scheme, and 
to identify in the latter formalism what are the different baryonic sector 
contributions to the partition function, the simplest way is to introduce a 
chemical potential \i since it couples directly with the quark number operator 
N. The Grand Canonical Partition Function (GCPF) is given then by the 
following expression 
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Z GC = Tr ( e -*C*-*")) . 



(4) 



Using now the previous decomposition of the Hilbert space in different 
baryonic sectors we can write the GCPF as 



which is the standard decomposition of the GCPF as a sum of canonical par- 
tition functions. Setting the chemical potential \x = we recover expression 



The standard way to introduce a chemical potential in lattice regularized 
QCD is to multiply all temporal link variables pointing forward (backward) 
by e M (e _M ) 0, [[|. The previous decomposition of the GCPF as a sum of 
canonical partition functions corresponds on the lattice to the Polyakov loop 
expansion of the partition function 



Z G c = a3V x e 3V *T + aa^-ije 3 ^* -1 ^ + + a + + a^ Vx e' Wx T , (6) 



where V x is the lattice spatial volume, the temperature T is the inverse lattice 
temporal extent and the indices in the sum are numbers of vanishing triality 
as a consequence of the fact that physical states are made up of a multiple 
of three number of quarks. This last property can also be derived from the 
invariance of the pure gauge action under Z(3) global transformations. 

The relation between the integrated coefficients in (6) and the contri- 
bution cik{U) from a single gauge configuration is 




(5) 



k 



(3). 



4 



a k = J[DU]a k {U)e'^ u \ (7) 

where Sq(U) is the standard pure gauge Wilson action, (3 the square of 
the inverse gauge coupling and a,k(U) is the contribution to the determinant 
of the Dirac-Kogut-Susskind operator containing a net number of Polyakov 
loops equal to k. 

Expression (6) for the partition function has been written assuming the 
Kogut-Susskind regularization for the fermionic degrees of freedom. Under 
such a regularization, 3V X is the maximum number of Polyakov loops which 
can appear in the fermion determinant for a lattice with V x spatial points. 
Equivalently 3V X is also the maximum number of three-colored quarks of the 
Kogut-Susskind type we can put in a finite discrete space with V x points 
following the Fermi-Dirac statistics. 

From the comparison of expressions (5) and (6) we can conclude that the 
coefficient a k is equal to the canonical partition function Tr k (e~^) computed 
in the baryonic sector of the Hilbert space having a baryonic number equal 
to | 0. 

In order to analyze the relevance of the different coefficients a k for different 
values of the baryon number |, we will consider the fermion determinant for 
a given gauge configuration in presence of an imaginary chemical potential 

detA(77,C7) =5> fe ([/)e ifc ", (8) 
k 
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where the sum in (8) is over all integer numbers between —3V X and 3V X . 

The antihermiticity and 75 anticommuting properties of the massless 
Dirac operator are preserved after the inclusion of an imaginary chemical 
potential. It follows that the positivity of the fermion determinant is also 
preserved and 

det A(r), U)>0 (9) 

for any gauge configuration and any rj. This property of positivity will be 
fundamental in order to derive bounds for the partition functions dk which 
will allow us to get interesting conclusions on the relevance of the different 
baryonic sectors. In a similar way Vafa and Witten used the positivity of the 
effective fermionic action in vector-like theories to get bounds which drive 
to the impossibility to break spontaneously vector-like global symmetries as 
isospin or baryon number conservation ||. 

A simple Fourier transformation in (8) gives 

<*k{U) =7^-1 e- lhv det A{r], U)d V (10) 

and, making use of the positivity relation (9), the following inequalities can 
be derived 

ao(U)>0, \a k (U)\ < a (U) Vfc. (11) 

These relations, which hold for any gauge configuration, imply for the 
integrated coefficients or canonical partition functions dk the inequalities 
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The relation (|T2| ) between canonical partition functions induces interest- 
ing physical consequences in the thermodynamics of this model. To analyze 
them let us consider the free energy density, 



/ = £logZ. (13) 

'X 

Using the obvious relation = d_k the partition function at zero chem- 
ical potential can be written as 



a relation which gives for the free energy density the following expression 



The first contribution in fllBD is the free energy density of the vanishing 
baryon number sector. The second term contains the contribution of all the 
other non-vanishing baryon number sectors. Relation (12) implies that the 
only contribution surviving the thermodynamical limit is the first one. The 
free energy density of QCD at any physical temperature T is given therefore 





(14) 





(15) 



by 




(16) 
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The previous expression implies that all baryonic thermal fluctuations can 
be neglected in the thermodynamics of QCD (at zero chemical potential). 
Notwithstanding that, this result does not necessarily imply that the proba- 
bility pk to get a sector with non-vanishing baryon number vanishes. What 
happens is that if — = is a non-vanishing number in the infinite volume 
limit, sector k has the same free energy density as sector and therefore the 
same thermodynamical properties. 

The CPT theorem || states that the vacuum state is CPT invariant. 
Since the baryon number density is an order parameter for the CPT symme- 
try, the theorem enforces the probability distribution function of the baryon 
number density to be a 5 function centered at the origin. Our results are 
therefore consistent with the CPT theorem but also give some further infor- 
mation: the partition function of QCD at any physical temperature can be 
computed not only in the vanishing baryonic density sector but also in the 
vanishing baryon number sector i.e., all baryonic thermal fluctuations can be 
neglected. 

This is something different from the standard statistical mechanics result 
that relative fluctuations go to zero as the thermodynamical limit is ap- 
proached. In fact we state that all properties related to the free energy den- 
sity can be derived considering the zero baryon number sector only. Clearly 
the result does not apply to the derivatives respect to the chemical poten- 
tial i.e., the baryonic susceptibility at zero chemical potential can not be 
evaluated from the vanishing baryon number sector only. 
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2. Some interesting physical implications 

Once we have demonstrated that thermal fluctuations of physical states 
with non-vanishing baryon number can be neglected in the thermodynamics 
of QCD, we will devote the last part of this article to the discussion of some 
interesting physical implications of this result. 

i) The first interesting feature is the possibility of analyzing the deconfin- 
ing phase transition at finite temperature by means of an order parameter. 
The mean value of the Polyakov loop or Wilson line || is related to the en- 
ergy of an isolated quark in the fluctuating background at a given physical 
temperature T. The standard wisdom is that the mean value of this oper- 
ator is not an order parameter in the full theory with dynamical fermions 
since the determinant of the Dirac operator, which appears in the integration 
measure, breaks explicitly the Polyakov symmetry. The possibility of using 
the mean value of the Polyakov loop as an exact order parameter in the full 
theory with dynamical fermions, which was pointed out in early works [0], 
[§], becomes apparent in the formalism here developed. In fact a (U), the 
contribution to the fermion determinant containing a vanishing net number 
of Polyakov loops, is invariant under Polyakov transformations. Even more, 
ao(U) is also independent of the boundary conditions for the fermion field 
in the temporal direction. Periodic or antiperiodic boundary conditions give 
rise to the same thermodynamical properties. 

ii) Theoretical prejudices based on simple models tell us that the high 
density phase of QCD should be similar to the high temperature phase in 
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the sense that both phases would be described by a free quark-gluon plasma. 
However the physical background in the two limiting cases, high temperature 
at vanishing chemical potential and high baryonic density at vanishing tem- 
perature, shows important differences. In the first case the ground state is 
dominated by thermal fluctuations of physical states with vanishing baryon 
number whereas in the second case, physical states with a high baryonic 
density are the only relevant contributions. Therefore even if the theoretical 
scenario of a quark-gluon plasma phase describes actually the two phases, 
we should expect important phenomenological differences between the two 
limiting extremes of the (/i, T) phase diagram. 

iii) The results here obtained allow also to solve the inconsistencies found 
in quenched QCD simulations at finite temperature. As it was shown a few 
time ago ||, |T0] quenched simulations of QCD at finite temperature and in 
the broken deconfined phase gave inconsistent results for the chiral conden- 
sate and different values for the critical temperature of the chiral transition 
were found, depending on the phase of the Polyakov loop in which the dynam- 
ics settles. These results, which put in question the validity of the quenched 
approximation, are not unexpected. In fact the integration measure is in- 
variant under Polyakov transformations in this approximation but the chiral 
condensate operator is not. A value of the chiral condensate depending on 
the Z(3) vacuum in which the dynamics settles in the broken phase is a rather 
expected result on physical grounds. However, as previously shown, the free 
energy density of QCD can be computed neglecting all the contributions to 
the partition function of non-vanishing baryon number sectors. We can take 
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ao(U) instead of det A(Z7) in the integration measure of full QCD and con- 
sistently the logarithmic mass derivative of ao(U) instead of the trace of the 
inverse Dirac operator as chiral condensate operator. Within this prescrip- 
tion, the chiral order parameter is invariant under Polyakov transformations 
and its mean value is therefore the same in each of the three possible Z(3) 
vacuums of the quenched broken phase. 

iv) Another topic in which the results here reported could give rise to an 
interesting improvement is the one concerning finite size effects and Monte- 
Carlo convergence in numerical simulations of finite temperature QCD. As 
previously discussed, most of the contributions of the different baryon num- 
ber sectors, specially those corresponding to a non-vanishing baryon number 
density, will have a vanishing probability in the thermodynamical limit. They 
are pure finite size effects and a natural expectation is that with our prescrip- 
tion of taking ao(U) instead of det A(C7) in the integration measure, we will 
approach the thermodynamical limit faster. But even if some of these contri- 
butions have a non- vanishing probability, our prescription of taking ao(U) in 
the integration measure should improve the Monte-Carlo convergence in nu- 
merical simulations of finite temperature QCD. In fact the coefficient ao(U) 
has more symmetries than the determinant of the full Dirac operator and, as 
we have learned from the GCPF computation in the investigations of QCD 
at finite density JTTJ] , a (U) is less fluctuating than det A(Z7). This last 
feature is the manifestation of what is known as the sign problem in finite 
density QCD JD|, a problem that has delayed enormously the progress in 
this field. 
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The central coefficient ao(U) is positive definite for any gauge configura- 
tion but unfortunately this is not the case for almost all the other coefficients. 
Most of them fluctuate violently with a rather constant absolute value but 
with an almost randomly distributed phase. Things here are very similar to 
what happens in the unpleasant situation we find when try to compute a small 
number by averaging numbers which are several orders of magnitude larger 
than their mean value. We would need a prohibitive statistics, increasing 



exponentially with the lattice volume in the GCPF computations |L2|, |L3 



to get significant numbers. The contribution of most of the non-vanishing 
baryon number sectors to the partition function of finite temperature QCD 
at actual statistics behaves therefore like a noise, and a natural expectation 
is that convergence of any Monte-Carlo procedure should improve by killing 
this noise. 

The exact computation of a (U) is very hard. The best way at present is 
through the construction of the quark propagator matrix |TJ| and requires the 
exact diagonalization of a non hermitian matrix for each gauge configuration. 
However, as follows from the analysis here developed, there is a simpler 
way to implement things in practical standard simulations of QCD at finite 
temperature. It consists in the introduction of an extra abelian degree of 
freedom e tv coupled to all temporal links and having a dynamical character. 
The integration over this new dynamical degree of freedom will kill all the 
undesidered contributions to the partition function. 
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3. Wilson fermions 

We have shown in the previous sections of this paper how the thermody- 
namics of QCD, when regularized in a space-time lattice and using staggered 
fermions, is controlled by the contribution to the fermion determinant with 
no net number of Polyakov loops, i.e., by the thermal fluctuations of physical 
states with vanishing baryon number. 

The two main ingredients to get this result are the conservation of bary- 
onic charge in QCD and the inequalities (12) which tell us that the partition 
function at fixed baryon number reaches its maximum value in the Hilbert 
subspace corresponding to zero baryon number. The first one of the two in- 
gredients is independent of the lattice regularization for the Dirac operator. 
The second one however is based on the positivity of the determinant of the 
Dirac operator for any gauge configuration and any value of the extra-abelian 
degree of freedom e tr >. 

Since we have made use of the hermiticity and chiral properties of the 
Dirac-Kogut-Susskind operator, it is natural to ask whether our result applies 
to any fermion regularization or rather it is related to the presence-absence 
of the chiral anomaly. 

Let us say from the beginning that even if we have not yet a definite 
answer to this question, there are strong indications suggesting that the chiral 
anomaly does not play any relevant role here. These indications come from 
the analysis of the properties of the fermion determinant for Wilson fermions. 
As well known, the Dirac- Wilson operator A can be written as 
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A = I - kM, (17) 

where k is the hopping parameter and the matrix M verifies the following 
chiral relation 



75M75 = M+. (18) 

Equation (|18|) implies that if A is eigenvalue of M, A* is also eigenvalue 
of M, i.e., the fermion determinant is always real. However it could be 
negative and in fact this unpleasant situation has been found for some gauge 
configurations in numerical simulations of the Schwinger model done in the 
unphysical strong coupling region However the unitary character of 



the gauge group implies that all the eigenvalues are upper bounded by the 
relation 



|A|<8 (19) 

which implies that for n < 1/8, det A > 0. It is easy to verify that under the 
previous condition k < 1/8, the positivity of det A also holds in the presence 
of the extra-abelian degree of freedom previously introduced. 

In other words, all the results of this paper can be extended in a straight- 
forward way to Wilson fermions if we impose the restriction k < 1/8, i.e., 
the hopping parameter region associated to a positive bare fermion mass. 
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